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We discuss vibrational localization problems in glasses and disordered media in this paper. It is claimed that the essence 
of the localization problem is already observed in disordered lattice models. These kinds of vibrations belong to a different 
universality class than bonded electrons. Specifically, The eigenvectors are extended even in two dimensions. Moreover, the 
correlation exponent does not diverge in the transition from localized to delocalized states. Furthermore, the volume of the 
extended states is scaled according to the distance from the transition as V ^ \uj — loc\L'^ ■ Interestingly, boson peaks can be 
observed in the density of states in both two and three dimensions. We studied the eigenstates of this problem and analyzed 
the scattering effects in these lattices. Importantly, we found that in two dimensions the boson peak, the localization edge, 
and the beginning of anomalous scattering are at the same frequency. In three dimensions, however, there are three separate 
regions: (1) localized, (2) weakly scattered, and (3) anomalously scattered. Finally, we discuss the relevance of this study to 
actual experiments and glasses. 

PACS numbers: 63.20 Localized Modes, 63.50 Vibrational States in Disordered Systems, 73.23 Theories and 
Modes for Localized States 



I. INTRODUCTION 

It is well-known that a disorder in the properties of 
a physical system can in turn induce localization and 
scattering effects. For the past three decades, disordered 
electronic systems have been the main focus of much re- 
search, having received enormous attention as manifested 
by the numerous books and reviews on the subject [1-10]. 
In contrast, though the disordered vibration problem is 
as important as the electronic one, we have found almost 
no discussion on the analogous problem of phonon local- 
ization. Though the vibrational properties of disordered 
solids were reviewed in several papers |]lT| , subjects like 
localization, scattering, and properties of the eigenstates 
were not analyzed. There are various reasons for this. 
The main one is that localization effects for electrons 
are more pronounced physically for statistical reasons. 
They occur at a well-defined Fermi energy state, so it is 
much easier to observe them under experimental condi- 
tions. In contrast, the localization effects of phonons are 
often obscured behind the calculations and the statis- 
tics; therefore, their role in influencing properties like 
heat conductivity or scattering is much less clear. We 
believe therefore that the theoretical understanding of 
these localization effects is crucial to the understanding 
of glasses. 

Unfortunately, we found almost no systematic treat- 
ment of vibration in disordered media above one- 
dimension in the literature. However, there is some liter- 
ature relating to two main issues. The first is a bonded 
network with cuts in the connections, which is the perco- 
lation problem. In such a problem anomalous fractal ex- 
citations called fractons may appear instead of the usual 
long-range phonons. This phenomenon was first men- 
tioned in [p^ and has been repeatedly observed in ex- 
periments and simulations. The second type of problem 



comprises simulations of 'real' glasses Though such 
simulations do yield some information about the relevant 
systems, it is rather difficult to obtain any concrete ideas 
because of various reasons. For example, simulations of 
glasses are completely uncontrolled in as much as the 
source of the disorder and its size are unknown, since the 
glass is defined by a dynamical relaxation algorithm. It 
is indeed possible to study these glasses and probe the 
degree of elastic and structural disorder; however, the 
disorder is still uncontrollable. Clearly, the effects of the 
disorder are manifested in many experimental measure- 
ments in glasses, but it is difficult to make a reasonable 
and rational interpretation of the data. 

Disorder in glasses and other materials can originate 
from various physical sources. Structural characteristics 
can change the topology of the material and thus affect 
the vibrations. Another major source of disorder is elas- 
tic variability in the local elastic constants. In simple 
glasses this is the main source of disorder, which is caused 
by the quenched pressure that remains in the glass after 
the quench from the liquid. Such quenched pressures are 
correlated with the existence of second-order gradients 
that generate large variations in the local elastic con- 
stants 1 16 Another important component for disor- 
der is the mass disorder that occurs when two or more 
kinds of atoms exist in the lattice. 

In glasses there is a set of anomalous thermodynamic 
properties that do not appear in crystals |^,|9|. Exam- 
ples of such behavior are anomalous heat conductivity 
and excess vibrations at low frequencies termed boson 
peak in the literature [ p^|Jl9[| . Scattering mechanisms are 
thought to mainly contribute to the depletion of the ther- 
mal transport in glasses. A crossover in the scattering of 
the glass is expected to be observed on a characteris- 
tic scale which is usually identified with the loffe-Regel 
limit, where the scattering length is of the same order 
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as the wave-length [^3|,|4|. In fact we observe two ba- 
sic scales in three dimensions as will be discussed later. 
Both of them do not fit the I-R demand but they are in 
the same range. In the second part of this paper we focus 
on scattering effects. Since getting experimental 15 1 and 
numerical information is relatively difhcult [ p^ ; it is eas- 
ier to use lattice models to obtain information about the 
vibrational properties. Indeed, we shall show here that 
many of the basic phenomena observed in glasses are also 
observed in lattice models. 

In the one-dimensional case, it is known that the vibra- 
tions are localized [ p0pl| , |2^ . However, the localization 
length in these systems is scaled as the inverse square of 
the frequency. So at low frequencies, we observe phonon- 
like modes with a localization length much larger than 
their wavelength, and with a low frequency spectral den- 
sity of order one. In a few theoretical papers on this 
subject 24|, it was contended that the situation at 
higher dimensions is analogous to that in electrons. We 
present here sufficient numerical evidence to show that 
this is not the case. 

The vibration problem generally consists of a vector 
translation problem. A simpler version of the problem is 
the problem of scalar vibrations where there is a scalar 
translation at any point. This is a realistic problem in 
one and two dimensions (e.g., a membrane). However, 
scalar models are, as we shall see, interesting enough the- 
oretically that we should start with a discussion of such 
models. 

There are in fact, there arc three non-equivalent phys- 
ical systems which, under certain constraints would give 
the same dynamical matrix: (1) the vibration problem 
(phonons) in disordered media, (2) free electrons with 
randomness in their transmission coefficients, and (3) the 
random diffusion problem. The dynamical structure of 
these problems is similar when the following constraints 
are fulfilled: (1) the transmission coefficients are sym- 
metric and positive, (2) the relevant dynamical variables 
are scalar, and (3) the masses are constant. 

Variants from these constraints, i.e. the vectorial case, 
the negative connectivity, and the variable masses, dis- 
play different results and will be discussed in other papers 
p5| . However, the time dependence, the statistics, and 
physical properties of these systems are very different. If 
any of the constraints are broken, we can expect to ob- 
serve a different behavior, as will be described later in 
detail. 

In this paper we extend our discussion of an earlier, 
shorter paper on the subject In section II we de- 

scribe the concepts of localization in electronic theory, 
and in general. In section HI we demonstrate that the dy- 
namical matrices of the three previously mentioned sys- 
tems are the same, when the parameter constraints that 
we mentioned are satisfied. In section IV we discuss the 
spectra and localization of the scalar models. We first 
analyze the one-dimensional case and then proceed to 



analyze numerically lattice problems that are analogous 
to the models studied in the electronic case[l-10]. By 
increasing the lattice size we can characterize the eigen- 
states. The participation ratios of the states are used to 
measure the degree of localization. A localization tran- 
sition between extended and localized states is observed 
both in two and three dimensions. Furthermore, in two 
dimensions there are three kinds of modes, namely, ex- 
tended, localized, and multi-fractal separated by mobility 
edges. 

However, the three-dimensional transition is different 
from the bonded electron transition[l-10]. Specifically, 
there is no divergence of the correlation length near the 
localization edge at Wc, and the volume of the extended 
states (defined by the participation ratio) converges to 
zero at the mobility edge. It is analogous to the per- 
colation phase transition above its threshold, where the 
global percolation cluster becomes a finite fraction of the 
size of the system. Such a transition implies zero elec- 
tronic conductivity near the mobility edge, which can also 
have a major effect on the heat conductivity of phonons 
in glasses, where there are well-known anomalies (for ex- 
ample, a plateau in heat conductivity in medium tem- 
peratures), see for example Ref. [|l9|. 

We have analyzed the eigenvectors of the two- and 
three-dimensional systems. Importantly, we found that 
in the two- and three-dimensional systems studied there 
are several scattering regimes. The first is at low fre- 
quencies in two and three dimensions where there is no 
interaction between the modes, and they are character- 
ized by well-defined wave vectors. The upper cutoff of 
this range depends on the system's scale. Above this 
cutoff the mapping between the state and the wave vec- 
tor fails because of scattering. However, an absolute 
value of the wave vector is still defined for each state. 
The inverse scattering width is defined by the variance 
around this mean value. This second range is called a 
'normal' scattering regime. We believe that it extends 
to zero frequency when the system is enlarged. Above 
this range beyond a point cj;, we find anomalous scatter- 
ing. Above these three regimes there are localized states 
where the wave length is hardly defined at all. For nor- 
mal scattering we find that the scattering width scales as 
r ~ ^ji 'i'io.i jj-^ dimensions and F ~ ^^^s.sio.s three 
dimensions within a relatively small region oi lu. In two 
dimensions there is an overlap between the appearance 
of anomalous scattering, the boson peak, and the local- 
ization edge i.e. lOc — LOb. However, in three dimensions 
there are two distinct scales. The first scale is the local- 
ization edge, and the second one is located at the boson 
peak and at the anomalous scattering edge Wf,, where a 
drop-off in the participation is observed. 
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II. LOCALIZATION AND ELECTRONS 

We shall review briefly presented some background 
on localization and its relation to the electronic case. 
The classical example for such phenomena is an electron 
in a disordered potential. An electron that obeys the 
Schrodinger equation with spatial disorder in the voltage 
is known to be localized for a strong enough disorder. We 
examine the equation 

zV-j = AV-j+F.Vj-, (1) 

where j is the position, A is the Laplacian operator on 
the relevant lattice, and the potential Vj is a random 
function. The variance of V affects the character of the 
eigenstates, the conductivity, and its dependence on the 
system size[l-10]. 

In one dimension and for any type of disorder all elec- 
tronic states are localized: the wave-functions decay at 
large distances as 

< i}{r)^{r + ro) exp(-|r - ro\/0, (2) 

where ^ is the localization length, and r and rg are any 
given lattice points. This leads to an exponential local- 
ization of the conductivity. 

In two dimensions the situation is marginal and de- 
pends to a large extent on the models. For any disorder 
in Vj the wave-functions form a multi-fractal set of the 
system. The multi-fractal exponent is weakly dependent 
on the energy. Above a certain critical disorder all elec- 
tronic states are localized. 

In three dimensions both extended and localized states 
exist. There is a critical energy Ec that defines the 
transition from localized to extended states. This tran- 
sition is accompanied by a diverging correlation length 
£,{E) [E — Ec)~'^ , where ^' is a positive number. In 
three dimensions, states are extended below or above a 
mobility edge and localized below or above it. The elec- 
tronic conductivity is expressed as a ^ [E — EcY 
For electronic models with magnetic fields there is also 
a localization transition in the two-dimensional case 
Classical waves (photons and waves) were proposed to 
behave similarly (see and Ref. therein). 

There are various criteria for localization cited in the 
literature, the first of which is the asymptotic behavior 
of localized states. It is usually described in terms of the 
envelope of the wave- function as stated in Eq. 2. 

It is often helpful to consider the moments of the wave- 
function values. For a particular normalized ipr{E), a 
generalized participation ratio is defined as 

l'^E)^l/iY,\A\''iE)), (3) 

r 

where the summation is done over all lattice sites, and q 
is an arbitrary real number. This measures the portion 



in space where the eigenvector is markedly different from 
zero. The dependence of on L can be written as 

1^1 ^ L^'-i\ (4) 

where D{q) is a fractal exponent. is the well-known 
participation ratio. The values of the Ds become con- 
stant in the limit of an infinite system. If the states 
are extended, D(2) = d. For multi-fractal states J^' is 
scaled as L^^'^^ with D{q) < {q — l)d. For a locahzed 
state, D{q >0) = and J"* - 

Another useful measure of localization is the transport 
coefficients of the system. In the electronic and thermal 
states the transport is defined by the Kubo formulas or by 
transport definitions [8]. A localized electronic problem 
has zero conductivity at zero temperature, whereas in 
multi-fractal cases it is believed that the conductivity 
converges to zero. 

III. MODELS 

A. Vibrations 

We begin this section with a discussion of the phonon 
model. If there is a particle system that has a stable 
equilibrium, it is possible to displace it by Ui from the 
stable particle positions. If the us are small enough there 
is an elastic expansion around the equilibrium positions. 
Accordingly, vibration Hamiltonians of the following kind 
can be written as 

^ = ^E"^''«)' + ^ E *"A^x<, (5) 

where uf is the a component of the translation in site i, 
and $ represents the interaction terms. 

We shall limit our discussion here to scalar elastic 
movements in an orthogonal direction to the state. (This 
will inhibit the interaction between the transversal and 
longitudinal modes, which have a role in vector models). 
So the discrete elastic Hamiltonian can be given as 

H ^^^m^iilif + ^^K^j{u^-Uj)^, (6) 
i ij 

where Ui represents elastic movements and Kij represents 
the elastic constants. 

There are two obvious constraints on the parameters. 
The first is the symmetry Kij = Kji. The second is that 
H must be non-negative for all realizations of u. This 
is a non-local demand for the elastic constants that is 
satisfied by a positive Kij . The single atom dynamics is 
defined by 

miili = E -^y ("j - "j")' (7) 

3 
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where rrii is the particle mass. 

Assuming that Ui = exp{iujt)ui{u!) we obtain an eigen- 
value problem represented by 

-miUj'^Ui = ^^Kij{ui -Uj). (8) 

j 

Note that the eigenvalue problem is generally not sym- 
metric because of the existence of the rrii terms. 

B. Diffusion and 'Free' Quantum Particles 

Free random diffusion has a similar dynamical form. 
We define particle densities pi and random diffusivities 
Dij between the nearby sites. The currents are Jij — 
Dij{pi— pj) and we obtain the random diffusion equation 

Pi^^D.jip,- pj). (9) 
j 

Defining pi{t, p) = cxp(— /it)pi(/i) results in a similar 
problem with p = —uP'. Since negative densities in p are 
not allowed, should be greater than zero. In contrast 
to the phonon model, there is no condition for detailed 
balance, and D does not have to be symmetric. 

A free randomly moving quantum particle can be de- 
scribed by 

= ^M,j(?A, -Vj) (10) 

for arbitrary Hermitian Af^j. For real and non- negative 
A/s and for the linear solution = exp{iEt)-i[ji{E), we 
obtain a similar eigenvalue problem. 

To summarize, the random realizations of these three 
linear models have the same linear dynamical structure, 
but only under the following conditions: the transmis- 
sion coefficients Af, D, and K are symmetric, real, and 
positive, and the phonon model is scalar and possesses 
a single mass. All these models are Positive Random 
Laplacian models. However, the time dependence, statis- 
tics and other properties are different. In this paper we 
use, unless stated otherwise, the phonon terminology. 

The difference between these Random Laplacian mod- 
els and the electronic models is already manifested in one 
dimension |2^j2l|,^ . The single mass vibration model in 
one dimension can be mapped (see in the next section) 
to an electronic model with a disordered potential that 
decays with the square of the frequency. In higher dimen- 
sions, however, we are not aware of any solution. Since 
we know that this model is equivalent to the diffusion 
model, it is clear that for positive connection coefficients 
any density will diffuse to infinity. Therefore, the corre- 
lation length can either explode at low frequencies or the 
system can be extended. 



IV. DEFINITION AND SIMULATION OF THE 
MODELS 

A. Definitions 

Here we focus on an equal mass (m^ = 1) model in a 
cubic lattice, of d dimensions, with a unit distance and 
size L. We considered nearest neighbor interactions, with 
two elastic constants Ki and K2, with a ratio a = K2/K1 
between them, and a probability p to find Ki and 1 —p for 
K2. Since a < 1 is the same as a > 1, we will discuss only 
the former limit. The normalization Ki — 1 determines 
a scale of frequencies. This model was suggested in Q 
and simulated in |Q for small as. 

There are two simple limits for this model, namely, 
a — I, which is the ordered lattice limit, and a = 0, 
which is for percolation. In the ordered lattice limit, all 
states are extended and ordered. The percolation limit, 
however, is more involved. Above the percolation thresh- 
old, we normally find phonons on a scale much larger 
than the percolation correlation length and fractal 
excitations, named fractons are found below this length 
[]l2| . Below the percolation threshold, the system would 
be completely localized. For a > 0, a different behavior 
exists. With extremely small values of a, there are fr ac- 
ton states together with phonons |^^. This is not the 
case were the values of a are larger and p is not close to 
the percolation threshhold. 

We treated the eigenvalue problem numerically using 
the LAPACK package for real symmetric banded systems 
for the equal mass model. Specifically, we used a cubic 
system of size L, where the distance between particles 
is unity. Most of our simulations were performed under 
free boundary conditions. However, we verified that the 
simulations were not changed under periodic boundary 
conditions. For this propose, we computed a set of eigen- 
values uJiS, and eigenvectors denoted as Ur{uji), where r is 
the lattice position. To verify the accuracy, we checked 
some examples for which the answers were known and 
compared the results of other codes to our computation. 

B. The One-Dimensional Case 

The one-dimensional problem on an ordered lattice is 
a classic example . Eigenstates can be written as 

Uk — eyip{i{kr + uj{k)t)). (11) 

If the lattice is disordered the states can no longer be 
written as periodic modes and thus the problem becomes 
localized. This general problem was solved in |^ , ^ , p2| 
and it is localized with a localization length that is pro- 
portional to u!~'^. We thought that it would be useful to 
test our code on this system to check both the numerical 
code and the dependence of our results on the system 
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size. The solution of the theoretical model in one dimen- 
sion can provide some insight as to why this model is so 
different from the bonded electronic model; therefore we 
present it in detail in Appendix A. 

Considering the above, we simulated the one- 
dimensional problem. Figure 1 displays our results for 
the participation ratio, which is as noted previously, pro- 
portional to the correlation length. This curve also shows 
an Lo^"^ dependence of the correlation length. As the sys- 
tem size is enlarged, wider regions of such behavior can be 
observed in the curve. For a definite system size _L, there 
is a cutoff in the frequency that is defined by L ^ ^c"^- 
Below this cutoff the correlation length is larger than L. 
Since lOc ^ L~^'^ and the number of states is propor- 
tional to L, there are L"-^ free states in such a system. 
This can be observed in the plateaus on the right-hand 
side of the curve. 

C. Two-Dimensional Models 

We now analyze the previous model in two dimensions. 
As mentioned previously, for the electronic model in two 
dimensions, any degree of variation in V generates fractal 
exponents D{2). These effects are observed in system 
sizes between 16 and 50. A transition in scaling for low 
noise levels can be easily observed. 

In the two-dimensional positive Laplacian case we were 
unable to distinguish any failure of normal scaling for 
small values of noise. Therefore, unlike the electronic 
case, a small disorder has no effect whatsoever on the 
scaling described by the participation ratio. For smaller 
values of a, localization is indeed observed. In Figure 
2 presents numerical results of systems with parameters 
p = 0.4 and a = 0.1. Figure 2a shows the density of 
states divided by w. A boson peak with a center around 
0.6 — 0.65 is visible. Figure 2b presents as a function 
of Lu for different system sizes. At low frequencies there 
are extended modes whose participation ratio scales with 
the system size. There is a well-defined mobility gap at 
LOc, separating the extended and quasi- localized states. 
Below the mobility edge the participation ratio can be 
written as 

I'^iiu) - \u;-u;c\L'^. 

All the states near the mobility edge are extended, 
though their volumes may be very small. The density 
of states in the extended region is close to the classical 
density N{llj) ^ uj. To analyze these states, we calculated 
the exponent D{2). Figure 2c shows the dependence of 
D{2) on w. Five scaling regimes of D{2) are visible. The 
multi-fractality of these states can be verified either by 
calculating the dependences of I^' or calculating of the 
distribution function of Ur and from point-point corre- 
lations. The transitions in D{2) are sharp within our 



numerical accuracy. Thus, there is a sharp transition be- 
tween the multi-fractal and extended states in this model 
at ujc = 0.63 ± 0.03. Note that this is also the position of 
the boson peak. 

D. Three-Dimensional Models 

The same kind of model was also simulated in three di- 
mensions. Figure 3a presents the density of states divided 
by u!^. A boson peak is observed at w « 0.6. Extended 
states and a localization transition to a fully localized 
state were observed (see Fig. 3b). Because of insufficient 
statistics, it is unclear whether there are quasi-localized 
states for this problem. Again, the scaling of the ex- 
tended eigenstates near the localized regime behaves as 

\UJ ~ LUc\L^ . (12) 

The states in the extended regime behave similarly to 
those e in two dimensions. Note that in this model the 
states up to a scale of w « 0.65 have a flat global scaling 
(see in Fig 3b). 

The phase transition from localized to delocalized 
states is different from that of the electrons. Specifically, 
there is no divergence of the localization length, as is ob- 
served in the electronic models, so ^ ^ (cj — uJc)'^. This 
also indicates zero 'conductivity' below the transition and 
will have an important effect on low heat conductivity in 
glasses p^ . The phonon volume in the extended regime 
shows a second order behavior similar to that observed 
for the volume of the infinite cluster in percolation mod- 
els. 

Let us now compare in detail comparison of our numer- 
ical results with the theory. Note that both theoretical 
works [^,^ are based on approximations. Here, com- 
parison of a numerical simulation with theory is possible, 
even for a system with a limited scale. The basic con- 
cept used in our simulations is the following: if there is 
an increase in the system size, localization lengths that 
are not dependent on the system size will appear in the 
simulation when the system scale becomes large enough. 
This is indeed what was observed in our one-dimensional 
computation. As the system size is increased, there is a 
decreased lower zone of frequency where participation ra- 
tios depend on system size, and an increased range where 
the they fit the theoretical predictions. This effect does 
not occur in two and three dimensions. There is a clear 
localization transition in two dimensions, but the states 
above it, even those with the smallest normalized partic- 
ipation ratios, continue to scale with the system size. A 
two-dimensional scaling length ^ ~ exp(— w~^) was pre- 
dicted in This indicates that above a cutoff given 
by L ~ ^, all states would scale with the system size and 
there would be a shift in the cutoff when L is increased. 
Both effects were not observed in our calculations. 
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V. SCATTERING AND EIGENVECTOR 
STRUCTURE 

A. The Two-Dimensional System 

As observed before the scattering of phonons and the 
shape of the scattering is an extremely interesting issue. 
We can expect to observe at least three scaling ranges in 
the shape of the spectra. Figure 4 presents shapes of var- 
ious wave-functions (left) and the square of their Fourier 
transform r(fc, w) (right). Figure 4a shows for lo « 0.1 a 
well-defined periodic mode together with a well-defined 
peak in Fourier space. As the frequency is increased, 
there is an increase in scattering. This is manifested in 
the creation of local domains with different directions 
and the existence of rings in Fourier space. At higher 
frequencies (w « 2), the shapes become localized. Fig- 
ure 5 gives examples of the radially integrated Fourier 
transform F{\k\,ix}) for the shapes presented in Figure 4. 

To quantify these effects, we measured the peak values 
of fco: which is scaled linearly with lo below the localiza- 
tion edge. Therefore, as stated before, they appear to be 
phononic. Above this edge over a limited region we find 

that OJ - fc0.85±0.3^ js^g^i- 

we assumed that the Fourier 
transforms can be represented as 



(fc-fco)2 + r2- 



(13) 



We make a fit with this expression to try to determine 
the amount of scattering. 

Fits with this expression are reasonably good up to a 
value oi uj K 1. Figure 6a shows the linear dependence 
of kf) at a frequency up to ujh ~ 0.65. Above this range 
the Lorentzian fits are less good because of the existence 
of larger tails at large fcs. In this range we observe (see 
Fig. 6a) a transition to an anomalous dependence of the 
wave vector uj ^ /j0.8±o.05^ Figure 6b shows the width F 
as a function of uj. In the smaller range scale, for small 
frequencies up to a scale of a; = 0.2 ± 0.05, the width is 
constant, since the system is not large enough to create 
scattering at these wavelengths. Above this value up to 
Wf, = 0.65, we obtain a dependence of F with /i, 

where /z = 1.7±0.2. In the second frequency range there 
is a smaller exponent, /i = 1.1 ± 0.1. 

This transition approximately overlaps the change be- 
tween extended and multi- fractal modes : = Wc- For 
the higher frequencies, the states are more localized and 
fits to Eq. O are not good. 



B. The Three-Dimensional Case 

We performed the same type of fit in this case as in 
the two-dimensional case. The results are presented in 
Figure 7. Here the range in Fourier space is smaller (one 



decade) and this makes some of our estimates inaccurate. 
As in the two-dimensional case, we observe for small uj 
values a range of constant peak width. For larger values 
of UJ there is a range with a very sharp rise, which can 
be fitted to F ~ (w - uj^f-^^°-^. For cjf, - 0.66 there is 
a crossover to anomalous scattering. When uj is further 
increased, there is a transition to localized states. We 
noted that this usually occurs in glassy models where 
the transition to a strong scattering transition happens 
within the extended range. 

Note that the dependence on F that we would obtain 
in optical Rayleigh scattering is F ^ cl'''+^, though this is 
not the same type of problem for the phonons ||15| . Our 
results in two dimensions are not consistent with this. 
In three dimensions they might be consistent with such 
scattering. 

In two dimensions the transition between normal and 
anomalous scattering is exactly the point where the lo- 
calization transition and the boson peak is observed. All 
these effects are clumped together in two dimensions. In 
three dimensions there are two transitions: the first is the 
localization edge ujc and the second is at the transition 
between normal and anomalous scattering and the boson 
peak at Wf,. This is a strong indication that in three di- 
mensions, the strong scattering effects begin before the 
localization edge effects do. This is also in accord with 
the commonly accepted situation as indicated by the lit- 
erature on glasses. 



C. Relation to Experiments and General Discussions 
on Glasses 

In actual experiments the dynamical correlations are 
measured by 

/oo ^ 
dt e-xjp{-iujt) < ^exp(zfc- {n{t) -rj{0))) >, 
-oo „■ 



(14) 



where the average is done in time at the relevant mea- 
surement temperatures. If we have a complete knowledge 
of the eigenvectors we can transform this expression to 

5(fc,c^) = 4,^<5(c^-c^0^('^,7^) 
I 

J2 e[ ■ ke] -keMiki^ " ^)) exp(-Wij) 



ij 

+h.o.t, 



(15) 



where the es are eigenvectors, N{oj, T) is the distribution 
function and the Ws are the Debye- Waller factors. The 
basic components in this sum are 



Gik,LJ,) = k^F{k,LJ,). 



(16) 
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Usually the sum is done over some range in w, so it repre- 
sents some average over the modes. In the normal scaling 
region the results we get from both kinds of fits are in 
accord. In two dimensions the fits can fail because of the 
existence of a longer tail in the Fourier transforms. 

In this paper we showed that indeed as speculated 
in the literature there in 3 dimensions there is a cross 
over frequency between anomalously scattered waves and 
weak scattering in which the width of the square of the 
fourier scales as (jj3-5±o.5^ rpj^^ scale of this crossover is 
not the loffe-Rogel scale so we don't know if it makes 
any sense to use this terminology. The conclusions from 
such a theory to transport is not clear. The reader can 
read some thoughts about it and some refrences from the 
glass literature in [?]. The transprt situation for states 
below the cross over should be discussed using different 
theoretical measures then ours. Nonlinear effects cannot 
be approximated using our method. 

We note that if there are nonlinear scattering effects 
associated with vibrations at low temperatures (like two- 
level systems) , they will be observed in scattering experi- 
ments but not in harmonic analysis. Furthermore if there 
are lifetime effects they will not be observed in such an 
experiment. 

As noted before, there is no actual available data for 
the normal scaling range, neither in real glasses nor in the 
numerical simulations, because of various experimental 
and numerical difficulties, so it is impossible to compare 
our results with experimental data. We hope that in due 
time this will be possible. 

VI. CONCLUSIONS 

In this paper we focused on the relatively simple case 
of the scalar vibration model. This can be demonstrated 
experimentally, for example, on a membrane. There are 
nontrivial predictions on the localization in such a sys- 
tem. We observed that the eigenvectors for these models 
are strongly scattered and that there are strong scatter- 
ing regions in the extended regime when the problem is in 
two and three dimensions. The dependences of the width 
of the Fourier transforms of the modes do not seem to 
be in accord with naive Rayleigh scattering estimates in 
two dimensions. 

We find that in the two-dimensional system there is 
only one critical scale that is observed by a localization 
edge, a boson peak, and a shift in the scattering and 
the dependence of uj on k. On the other hand, in the 
three dimensional system there are at least two scatter- 
ing scales and three scattering regimes. The first one 
is the localization edge and the second is the scattering 
edge where scattering becomes large. The second scale 
overlaps with the boson peak. A similar picture about 
Silicon glass is given in even though we do not agree 



with many of the statements in this paper about the na- 
ture of the states. This subject will be discussed in more 
detail in another paper where the transport coefficients 
will be calculated. 

Introducing of vector disorder and structural disorder 
will induce further complications, consisting of interac- 
tions between longitudinal and transversal modes. The 
existence of mass disorder is another strong mechanism 
for introducing more scattering, especially in the optical 
bands This subject is especially interesting because the 
representation of systems with mass disorder is equiva- 
lent to a Schrodinger equation with a correlated poten- 
tial. This will be discussed in another paper [^. An- 
other subject that remains largely untreated in our paper 
is the definition of electrical and thermal conductivities. 
Kubo and transport formulas can be used to estimate the 
dependence of these conductivities on temperature and 
elastic disorder this will be done in another study. 

We demonstrated that a set of models has a different 
behavior with the effect of noise than bonded electrons. 
This is manifested in two main ways: (1) the phonon 
models seem to be delocalized in two dimensions, and 
(2) they have a percolative behavior as a function of the 
frequency. The scale of the mode diverges as the distance 
from the critical points times the system volume. 

We focused on the relatively simple case where the 
movements are scalar. This can be realized experimen- 
tally, for example, on a membrane. There are nontrivial 
predictions for such a system. Introducing of vector dis- 
order and structural disorder will induce further compli- 
cations, consisting of interactions between longitudinal 
and transversal modes and possibly soft modes (to be 
discussed elsewhere). The existence of mass disorder is 
another strong mechanism for introducing more scatter- 
ing, especially in the optical bands. Notably, the eigen- 
vectors for these models are strongly scattered and there 
are strong scattering regions in the extended region in 
two and three dimensions. 

Another subject that remains largely untreated in this 
paper is the definition of electrical and thermal conduc- 
tivities. Kubo formulas can be used to estimate the 
dependence of these conductivities on temperature and 
elastic disorder , which will be done in another study. 

We acknowledge discussions with 

O. Gat, J.P.Eckmann, S.R. EUiott E. Courtens and R. 
Vacher. Ill for discussions of the subject. 
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APPENDIX A: THE ONE-DIMENSIONAL 
VIBRATION PROBLEM 

It was shown in |^,^ that the phonon model is equiv- 
alent to the electronic model with a rescaled interaction 
ViUJ^^. We shall briefly repeat the arguments because, 
first, we can use them to solve this problem and second, 
because it is reasonable to assume that there is an anal- 
ogous effect at higher dimensions. 

The one-dimensional equation can be written as 



Up'Ui 



) + A'i_i(ui_i - Ui). 



(Al) 



This model can be solved through the following transfor- 
mations: if we define Qn — Ki(ui — Wi+i), the equation 
is transformed to 



Q 



n+l 



(A2) 



where 51 



This is a Schrodinger equation with a 



potential that scales down with the frequency. Eigenval- 
ues explode exponentially as 



Qn = exp(n7(ri)). 



(A3) 



An equivalent equation for i?„ 
ten as 



Qn+i/Qn can be writ- 



R„ 



2 + niKi - 1/R„ 



(A4) 



We can now write an equation for the probability distri- 
bution P{R) in the steady state. Here we get 



The localization exponents are defined as 



7(f7) = j P{R)ln{R)dR. 



(A5) 



(A6) 



For small values of couplings or frequency the solution 
is defined by a fixed-point solution. Using the average 
value of the map h = fiA, {A = {p/Ki + (1 -p)/K2)) we 
get 



R=l 

which can be expanded to 



i?- 1 = 1 + 6^ + - + 0(6^/2). 
The localization length follows 



(A7) 



(AS) 



(A9) 



These results can be tested in a simulation as we here, 
mainly, to test the simulation. 
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FIG. 1. Log-log scale of the participation ratio (proportional to ^) as a function of system size for one-dimensional systems 
with a = 0.1 and p = 0.4, for sizes 128-8192 with jumps of factor two. We observed that the curves decay as ui~'^, as the bold 
line illustrates. As the system size is increased, this scaling range also increases. Note that the upper cutoffs are defined by 
L ~ ujc'^. Below ojc, is proportional to the system size. 
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FIG. 2. (a) The density of states A'^, for p = 0.4 and a = 0.1, divided by w for a system size L = 100. The peak is the extra 
number of states relative to the Uncar normal distribution, which is called the boson peak, (b) The participation ratio is /■*. 
The curves are averaged over equal intervals of ui of size 0.05 to reduce the /■* fluctuations. The extended modes are seen below 
uj ~ 0.63. (c) The scaling exponent D{2) is shown as a function of the frequency u!. The curve is averaged over intervals of 
frequencies lj in order to reduce the noise in 
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FIG. 3. (a) The density of states A'^, for p = 0.8 and a = 0.1, divided by uj^ for a thrcc-diniousioual system with a size of 
L = 22. The peak is the extra number of states relative to the linear normal distribution, which is called the boson peak, (b) 
The participation ratio normalized by with the same parameters. The solid line is for L = 8, dotted for L = 16, and 
dashed for L = 22. The curves are averaged to reduce the fluctuations. For the system with L = 8, we ensemble the average 
of 6 runs to reduce fluctuations. The extended modes are observed below w ~ 1 by their collapse into one curve. 
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FIG. 4. We present the two-dimensional shapes and Fourier transforms of the states of the 100'^ system with p = 0.4 and 
a = 0.1 for a; = 0.1(a), 0.2(6), 0.6(c), 1.0(d), 2.0(e). 
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FIG. 5. The radially averaged squares of the two-dimensional Fourier transforms F{k) for the previous uis are shown. 
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